Positive P-Representations of the Thermal Operator 
from Quantum Control Theory 



John A. SidlefE 

Quantum System Engineering Group 
University of Washington, School of Medicine 
Box 356500, Seattle, WA 98195, USA 
(Dated: January 26, 2004) 

A positive P-representation for the spin-j thermal density matrix is given in closed form. The 
representation is constructed by regarding the wave function as the internal state of a closed-loop 
control system. A continuous interferometric measurement process is proved to einselect coherent 
states, and feedback control is proved to be equivalent to a thermal reservoir. Ito equations are 
derived, and the P-representation is obtained from a Fokker-Planck equation. Langevin equations 
are derived, and the force noise is shown to be the Hilbert transform of the measurement noise. 
The formalism is applied to magnetic resonance force microscopy (MRFM) and gravity wave (GW) 
interferometry. Some unsolved problems relating to drift and diffusion on Hilbert spaces are noted. 
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Given spin-j operators s ~ (si, S2, S3), it is well known 
P, 0,13, 3 that a coherent state \x) can be associated with 
each unit vector x, such that {x\x) = 1, (£|s|&) = jx, 
and a resolution of the identity operator I is 



I = 



2j + l 
47r 



Sx \x) {x\ 



(1) 



47r 



Given an operator pj on each member of a sequence 
of spin-j Hilbert spaces j G {0, 5, 1, • ■ •}, we define a P- 
sequence to be a set of functions {Pj{x)}, and similarly 
a Q-sequence to be a set of functions {Qj{x)}, satisfying 



{x\pj\x) 



2j + l 
47r 

0,(A), 



(fxPj{x)\x){x\, (2a) 



4ir 



(2b) 



The properties of P-sequences and Q-sequences will be 
the main topic of this article. 

Q-sequences pose few mathematical mysteries because 
(|2b(l suffices to evaluate Qj{x) directly. In contrast, lit- 
tle is known about P^-representations, particularly about 
their positivity properties. The main obstruction is that 
Pj is nonunique due to the identity 



/ (fx Y^{x) \x){x\ =0 for Z > 2j, 

J47r 



(3) 



where Y^{x) is a spherical harmonic. To circumvent this 
P-representation ambiguity we will apply methods from 
quantum measurement and control theory. 

We will study the P-sequence of the thermal operator. 



„th 



P 



exp(— /3t's) , 



(4) 



where t is the thermal axis. The Q-sequence of p**^ can 
be computed by direct evaluation of l)2b|l. and the well- 
known result 111 l3| is 



Qj{x) = (cosh — i-t sinh 5/?)^'' 



(5) 



Our main result is: 

Theorem 1. A positive P-sequence for the thermal op- 
erator is Pj{x) — Cj/Qj^i{—x). 

The normalization coefhcient Cj is readily found from 
(|2a|l by taking the trace, with the result 



2sinh((j + i)/3) 



(6a) 



We first prove Theorem ^ nonconstructively by treat- 
ing it as an ansatz. We take the thermal axis t = (0, 0, 1) 
such that {j,m\p^^\j,m') = 5mm'&~^"^ ■ Taking matrix 
elements of (|2a|l . Theorem Q] is proved if Pj satisfies 



-/3m s 



2j + l 
47r 



ffxP,{x){j,m\x){x\j,m'). (7) 



47r 



We write \x) = D{(l),9,0)\j,j), where the rotation 



(8) 



carries t — (0, 0, 1) into x = (sin 61 cos (p, sin 6* sin 0, cos 9). 
By an identity due to Wigner 0,01 

(j,m|A) = (j,m|i?(0,0,O)b-,j) =i?^,/0,0,O) 

= 0-m)^^^e-*™'^cos^+™i6l sinJ-^i^, (9) 

we transform (O into an integral representation of the 
hypergeomctric function 2 Pi. Then via the identity 

2Pi(2 + 2j,l+j + m,2 + 2j,z) = (1 - z)-(^+™+i). (10) 

Theorem ^ follows immediately. 

But where did the ansatz come from? Why do positive 
P-representations even exist for thermal operators? 

We will now give a constructive proof that answers 
these questions. Our notation and methods are adapted 



2 



from Handbook of Stochastic Processes ■ Purely alge- 
braic details are not given, but when non-obvious idioms 
or strategies are employed, we will show them. 

The proof consists of three physically-motivated steps. 
In the first step, we will regard {ip) as the internal state 
of an open-loop dynamical system, and we will specify a 
sensor that monitors the spin axis. 

In the second step, we will install a controller that 
uses the (imperfect) sensor measurements to align the 
spin axis with the thermal axis t. We will prove that the 
control noise is precisely equivalent to a thermal bath. 

In the third and final step, we will construct Ito and 
Fokker-Planck equations for the observed states. Theo- 
rem lUthen emerges constructively, with Pj as the solu- 
tion of a Fokker-Planck equation. 

To begin Step 1, we define an open-loop uniaxial spin- 
ometer with generator s to be the Markov chain of quan- 
tum states (see, e.g., 0) defined by 



^l^n)/^ 



"Pa; 



Pa 
Pb 



(V'niBtSlV'n) 



(11) 



with increment operators A and B given by 

A = [cos{es) + sin(6ls)]/V2 , (12a) 
B = [cos(6ls) - sin(6's)]/V2 . (12b) 

The generator s can be any Hermitian operator. 

Here Pa + ^'b = 1 is guaranteed by A^^ A + B'^ B ^ I. 
The measurement strength is governed by 0, and we will 
assume <C 1. Terms of order 0'^ and higher will turn 
out to be nonleading and we will ignore them. 

We define the operator variance A(s) to be 



It is readily shown that in general 



(13) 



. , , I = if |'0n) is an eigenstate of s, 
A„{s){ ' ^ (14) 

I > otherwise. 

For E an ensemble average, the variance increment 

5£;„[A(s)] = £;[A„+i(s)] - E\I\,,{s)\ (15) 

is found by straightforward stochastic analysis to be 

5£;„[A(s)]--402^[A2(s)]. (16) 

The sequence {i?[Ai(s)], i?[A2(s)], . . .} is positive by 
()14|l and decreasing by (|16|l . therefore it has a limit. 
This implies a vanishing increment (|16|l . which implies 
lim„_oo £'„[A2(s)] = 0, which implies lim„^oo A„(s) = 
for every Markov chain in the ensemble (except a set 
of measure zero). This proves 

Theorem 2. Open-loop uniaxial spinometers asymptot- 
ically einselect eigenstates of the generator. 



Here we have used einselect in Zurek's sense [3: 

Einselection [excludes] all but a small set 
[of states] from within a much larger Hilbert 
space. Einselected states are distinguished by 
their stability in spite of the monitoring envi- 
ronment. 

Now we operate three spinometers simultaneously, with 
generators (51,52,^3). We call this an open-loop triaxial 
spinometer. We define a spin covariance cr„ 

{o'n)kl = (V'nISfcSilV'n) " (V'n I Sfc | V'n) (V"™ | Si | V'n) (17) 

and similar to (|13|l and H15|l we define an n'th stochastic 
increment for the i'th spinometer 



(18) 



By straightforward stochastic analysis similar to 1)16(1 , the 
net covariance increment is calculated to be 

3 

ixbE,\(T\ = tr SE,^„[(t] = -46^ tri;[cr„.cr*] . (19) 

i=l 

Physical intuition then suggests that similar to H14|l 
if and only if \ipn) is coherent, 



tl (T„ -(T* 



> otherwise. 



(20) 



To prove this we define a spin vector x = (iplslip) / j , 
such that aim — cf-mi = i^imnjxn- Without loss of gen- 
erality we choose a frame in which x = (0,0, 3:3). We 
define the symmetric tensor aim = \ {<^im + <^mi), and we 
temporarily regard xj, and {akV, fc > as a set of seven 
arbitrary real numbers. Then solely by reordering terms 
we construct two algebraic equalities 



(21a) 



trcr-cr* = Pa + 2p\, + \pc+ ipd + ^pl, + \pc + \pf 
tra-a ^ \f +p^ + 2p^ + \p^+ jpA + \pl + \pi (21b) 

whose individual terms are 



2\2 



Pa = (CT33 - X^) 



^23 



Pi 



] -cr33 



Pc = (o-ll " 0-22)^ 

Pd = (tro- - CT33)^ - + 1)- asaf . 

By construction, each term is nonnegative for all [ip) and 
zero for coherent \Tp). In particular, p^ and pf vanish if 
and only if CT33 = {ip\sl\^p) = f and jxs = {ip\s3\tp) = 
±j, i.e., if and only if \i[) is coherent. Then (|20|l follows 
immediately from the nonnegativity of H21a|) . Similarly 
()21bp implies a related inequality for &: 



trcr-cr ■ 



= if and only if \i[i) is coherent, 

> ^j^ otherwise. 



(22) 



Then by the same reasoning as for Theorem [21 1(1 9|l 



implies that lim„ 



, tr (T„'cr* — for every chain in the 



ensemble (except a set of measure zero). This proves 
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Theorem 3. Open-loop triaxial spinometers asymptoti- 
cally einselect coherent states. 

This completes Step 1 of our program. 

Remark: it seems reasonable that all Lie groups might 
asymptotically einselect coherent states |3, but the au- 
thor has studied only SU (2) P-representations. 

Now we begin Step 2, and focus on control and ther- 
modynamics. For t the thermal axis defined in we 
modify the spinometer matrices (|12b|l such that 

Al = e-*"(*x^)'» [cos(e'sfe) + sin(6lsfe)]/ V2 , (23a) 
Bl = e+'"(*><^)'« [cos(6lsfe) - sm{esk)]/V2 . (23b) 

We will call this a closed-loop triaxial spinometer with 
unitary feedback, because the operators exp(±?at x s) 
act cumulatively to align the spin axis with t. 

Closing the control loop does not alter the coherent 
einselection because the sole effect of a post hoc unitary 
operator on cr„ is a spatial rotation. Since tr cr„'(T* is a 
scalar, (|19f) still holds. Thus we have 

Lemma 1. Closed-loop triaxial spinometers with unitary 
feedback asymptotically einselect coherent states. 

An ensemble of closed-loop spinometers has a density 
matrix sequence {pi, p2, ■ • •} whose increment is 

3 

Spn - (^fcVn^fe + Bl^p,,Bl - p„) , (24) 
fe=l 

and we readily compute that 6pn — for p„ = p'^, with 
the thermal operator defined in |0J, provided 

a = - tanh i/3 or 1/a = - tanh jf3 . (25) 

We will show later on that p**^ solves Spn = uniquely, 
because the Fokker-Planck equation for p has a unique 
stationary solution. This proves 

Theorem 4. The density matrix of an ensemble of 
closed-loop triaxial spinometers with unitary feedback is 
asymptotically thermal. 

To connect (|24() with the thermodynamic literature, 
we set t — (0,0,1) and expand to order 6'^. The re- 
sult is equivalent to a thermal model given by Perelomov 
(per eq. 23.2.1 of Q). Gardiner gives similar model, (per 
eq. 10.4.2 of 0)- In Lindblad form we find 

Spn = -\l{v + 1) (s+s_p - 2s_ps+ + s+s_) 
— \"fu [s^s+p — 2s+ps_ + s_s+) 
+ 0'^ {S3S3P - 2s3ps3 + S3S3) , (26) 

where s+ = (si + is2)/V^ and s_ = (si — is2)/V^, and 
we have adopted Perelomov's variables 7 = —4a^9^ and 
v = — (l-|-Q;)^/4a. This completes Step 2 of our program. 



Now we turn our attention to Step 3, and focus on Ito 
and Fokker-Planck equations. 

The following idioms lead quickly to Theorem ^ We 
define a data three-vector d„ = {di ^, d2.n1 c's.n) by 

^ _ 1+1 for IV-n+l) CX IV'n) , 

''"~\-l for cx S,^ |^„) . 

Thus {di, d2, . . .} is binary- valued data. The mean 

E[dr.] = (V-n l^^fcU^fc I Vn) " | V'n) (27) 

satisfies 

E[dn] - 20E[Xn] = gsE[xn] (28) 

which defines gs = 20 as the sensor gain. We remark that 
dji/gs is therefore an unbiased measure of x„. 

We next define a zero-mean stochastic variable Wn by 

dn 5s Xn + Wn , (29) 

such that (to leading order in 0) Wn has the second-order 
stochastic properties of a discrete Wiener increment: 

E[iWn)k{Wn')k'] = Snn'Skk' ■ (30) 

An identity valid for a coherent state, 

{i>n\skSl\'lpn) = ^jSkl + j{j - \)xkXj + \ieklraXm , (31) 

gives rise to an Ito increment of conventional form 

5Xn = Xn+l - Xn = gla{Xn) + gs b{Xn)-Wn , (32) 

where a is called the drift vector and b is called the dif- 
fusion matrix. 

As a check, Lemma ^ requires that the stochastic mo- 
tion of Xn be confined to a unit sphere, and so does the 
coherent state identity (|31|l . since it is inhomogenous in 
X. For consistency, therefore, the increment of the m'th 
radial moment must vanish for all m when \x\ = 1. This 
increment is readily calculated have the general form 

(5£'„[|a;|'"] cx ^ni{m ~ 2)[a;„-b(a;„) • b''{xn)-Xn] , 

+ m\Xn\'^[^tl-b{Xn) ■ b''{Xn) + Xn-a{Xn)] . (33) 

Upon computing a{x) and b{x) explicitly we find 

a{x) = Ix [a{l - 2j)x-i- (1 + ^a^)] 

+ it [ a{l + 2j) - \o? x-i ] , (34a) 
b{x) = \[I - X ® x + a{i® X - x-il)], (34b) 

such that the radial increment 133|) indeed vanishes. 

A Fokker-Planck equation for Pj{x) is readily found 
from (|32|l and H34al -b') . Setting z = x-t we obtain 

= -^[a(l + z^) + 2ja(l - z^) - z{l + a^)]Pj{z) 

+ l&ii^ - - + ^')]Pji^) . (35) 
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which has the unique properly normaHzed solution 



1 - 



-, 2j+2 



1 — 2a x-t + 



(36) 



Substituting a = — tanh|;/3 per l|25|l yields Theorem ^ 
This completes the third and final step of our proof 

We now discuss some practical implications for exper- 
iments in magnetic resonance force microscopy (MRFM) 
and gravity wave (GW) interferometry. We first write 
the Ito equations H32|l in Langevin form by substituting 

bxn^ lldt'x{t') a{xn) rj*dt'a{x{t')) (37a) 
W,,^rg4^dt'x''{t') b{xn) ^ bix{t)) (37b) 

where St — 1/r is an interval and x^{t) is white noise 

E[x^{t)xl{t')]=Skk'S{t-t')/g!r. (38) 

Then taking d/dt, the resulting Langevin equation is 

X = rg^[a{x) + b{x)-{x^ ^ x)] , (39a) 

where x^{t) — x{t) + x^{t) is the measured spin axis. 

We see that x{t) is dynamically attracted toward the 
measured axis x^{t). Even open-loop spinometers ex- 
hibit this attraction, since for a = we find 

x\a=o ^ rg^l-jx + ^{I-x(E) x)-{x^ ~ x)] . (39b) 

Remark; a similar einselection-by-attraction is evident 
even in uniaxial spinometry, where it dynamically gener- 
ates the asymptotic "collapse" of \ipn) to an eigenstate, 
as described by H16|) and Theorem |21 

We now transform HSgal -bl to the second-order Newto- 
nian equation of an oscillator. To do this, we introduce 
a mass m and frequency ujq by defining 

q°P = (h/jmujo)^^^ [+S2 cos{ujot) - si sin(tJoi)] , (40a) 
p°P = {mujoh/jy/^ [~S2 sin(wot) - si cos{ujat)] . (40b) 

We verify that for states with z ^ 1 the canonical com- 
mutator [(7°P,p°P] = ifis^/j ~ ih holds. Defining the 
coherent oscillator coordinate q{t) to be 

q{t) = {j h/mLUoy^^{y{t) cos Uot — x{t) sin UJot) , (41) 

we find that (j39b|) takes the equivalent Newtonian form 



(42) 



where f°^^{t) is an arbitrary external force. From H38|l 
we obtain the noise spectral densities, which satisfy 



SqNfN{uj) — iiftsgncj. 
Here ^ JZodr e'^^^ E[q^ {t)f^ {t + r)]. 



(43a) 
(43b) 



We recognize ()43a|l as the standard quantum limit 
(SQL), but the expression ()43b|l for S^n^n is surprising. 
We will call it the Hilbert correlation because it asserts 
that f^{t) is proportional to the Hilbert transform of 
q^{t). Physically, it ensures that fluctuations in the ob- 
served position q^{t) are accompanied by force flucta- 
tions such that q{t) is attracted toward q^{t), as H39al -b) 
requires (and as is ubiquitious in spinometry). 

Is the Hilbert correlation observable? From H42|l it can 
be shown that it is undetectable in g^(t)'s response to 
the classical force f^^^it). Thus the standard quantum 
limit for classical force signals, as detected in experiments 
like GW interferometry, is not altered by the presence (or 
absence) of a Hilbert correlation. 

The author is presently investigating the conditions 
under which Hilbert correlations are observable; in spino- 
metric terms this seems to require a curved (nonlinear) 
dynamical manifold, as can arise, e.g., in MRFM when 
several spins are present. 

A major obstacle is that the differential geometry and 
topology of drift and diffusion functions on Hilbert spaces 
are poorly understood. The following would help: 

Problem 1. Given a finite- dimensional Hilbert space 
and on that space a set of Ito drift functions aQip)) and 
diffusion functions b{\%l))), either exhibit a set of spino- 
meter operators {Ai, Bi\ that generate a and b, or prove 
that no such operators exist. 

When students ask "How does the Stern-Gerlach effect 
work?" or "How does the standard quantum limit work?" 
the author increasingly answers them in terms of drift 
and diffusion functions. This focusses on the well-posed 
problem of designing these functions and realizing them 
in hardware, and yet makes clear how little is known, and 
how much remains to be discovered. 
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